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Some Infinite Sums Involving Zeros of J (x) 

Leo F. Epstein * 

(April 15, 1963) 
The limiting form for small </> values of the sums 

/v(4>)=£ (l/j8„r exp (-j8«*<« 

n=i 

_ " exp (-(3 n 2 (}>) 
M r «>) =£ <l/j8„) 2 '-> erfc (j8„ V*) 



A erfc W„V$) 

where r is an integer and /3„ is the Ait h zero of the Bessel function J (x), have been obtained. Related 
functions which arise in these problems are also considered. 

1. Introduction 

In the study of a physical problem involving diffusion processes in cylindrical geometry, a 
number of infinite sums have been encountered of the form 

FA<t>) - J (l//3„) 2r exp {-(i,?4>) (1-1A) 

n=\ 

GrW = J, /B„-y.0„) (1 " 1B) 

where /3 n is the nth zero of the Bessel function Jo(x). For these functions, the limiting behavior 
for small </> values is required. The technique of evaluation involves contour integration, the use 
of Laplace transforms, and several other methods of some interest, and the results given below 
are, except as noted, not known to be recorded elsewhere. In the cylindrical geometry problem 
in which these functions, and the related M r (4>) and Nr(<l>) introduced later, were found to occur, 
the <f> — > solutions examined correspond to the short time behavior of the system (see eq (4H3) 
below). It has been found that these quantities are frequently encountered in studies of diffusion 
and related phenomena, heat conduction for example. 
The asymptotic expression [1, p. 505] l for j3 n is 






Urn Pn^Pt^^Ctj 

^00 j = 



rr{n-\) 



1-2/ 



(1-2) 



*Vallecitos Atomic Laboratory, General Electric Company, Atomic Power Department, Pleasanton, Calif. This paper was first issued as a G. E. internal mem- 
naiulum, GEAP-3853, dated July 1, 1961. 

1 Eigures in brackets indicate the literature references at the end of this paper. 
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where the first few coefficients of the series are 

a =+l 

«!=+ 1/8 

a 2 = -31/384 

a 3 = + 3779/15360 (1-3) 

jmin is the value of the index for which the summand is a minimum, so that the series is terminated 
at the term just before the smallest [2, p. 140]. 

From (2), it is immediately apparent that the series in (1A) and (IB) converge, and that F r {(f)) 
and Gr(<f>) exist for all positive real values of r and 0. 

2. Evaluation of Fi(<£) 

To evaluate Fi((/>), consider the function 

C 

integrated about a contour 2 consisting of two lines parallel to the real axis and at an arbitrary 
distance, plus and minus, from it large enough to include z = ±ip v * . Applying Cauchy's theorem 
[3], 

^) =+ ^(Lim[^.^] + |nm[ 2(2 _^^.^]) (2 _ 2) 

where the two limits are the residues of the poles of the integral at zero and (3 n respectively. The 
first term here is iir\2 as can be found from the leading terms in the series for J n (z). By expanding 
the summand in a Taylor's series about z = /3 n , the second term can also be evaluated, and it is 
found that 

(2/i77)/(<M= 1 - 4 J (1/j© exp (- pl$) = 1 - 4F,(0). (2-3) 

Taking the Laplace transform, defined by [4] 

00. 

L P [g(<k)] = I g(<W exp {-p<}>)d4 

of both sides of (2-3) , 

L P [l-4F lW ] = (2MJ^L^^ 
c 

The integral here has a pole at z = ip* and the residue at this pole is readily found. The value of 
the integral, then, is 27ri times this residue [3], and 

LJLl - 4Fi(#] =+ 2p- 3 / 2 i-^-p- (2-4) 

where [5, p. 372] 

In(z) = i- n Jn(iz). (2-5) 

Using the asymptotic form [1, p. 203] for I n (z) for large z 

m . —1 

t m « e * Y ( _ 2z )-m nn + m+j) J_ 



2 The technique used here was suggested by Dr. G. M. Roe of the General Electric Research Laboratory, Schenectady, N.Y. (private communication June 13, 
1955). 
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it follows that as p — > °° 

L P I 1 - 4F,(0)] - + 2p-*l*-p- 2 - (l/4)p- 5 / 2 - (l/4)p- 3 + . . 

This form suggests letting 

l-4F 1 (<t>) = J 2a„4>"l 2 

n=\ 

and taking the Laplace transform of this expression [4, p. 137] 

L p [l-4F,(<p)] 

= (1/2) la ip - 3 l 2 + (1) \a- 2 p- 2 + (3/2) !a 3 p" 5/2 + (2) W" 3 + . . 



(2-7) 
(2-8) 



(2-9) 



Comparing the coefficients of like powers of p in (7) and (9) yields the values of the as and for 
small (f> 




(2-10) 



Note that, from this equation, using the notation of Watson [1, p. 502] 

(To |1)s 2 l/j6„ 2 = F,(0)=l/4 

n=l 

a well-known result. 

Differentiation of (2-10) with respect to (/> yields, for small </> 



(2-11) 







F (<f>) 






= (1/2) (1/ V^)</>-" 2 - 


-(1/4)- 


-(1/8) (1/ V^)</> ( 


1/2. 


-(1/16)0+. . . . 



(2-12) 



It is of interest to compare (12) with the similar expression [6, par. 6.491, p. 129] in which /?„ 
is replaced hy n 



J exp(-n 2 </>) = ( V^/2)</)-'/ 2 -(l/2)+ V¥</>-' /2 J exp (- rc 2 7T 2 /</>). 

n=\ n=\ 

The expression corresponding to (2-10) for small cj> is [7, pp. 234 and 484, and 19, p. 64]. 

J (l/rc 2 )exp(-n 2 (/>) 

H = l 

- + (7T 2 /6) - V^(/) ,/2 + ( 1 /2)(f} 



(2-13) 



-(<t>lirfl 2 2 (1M 2 ) [1 "(3/2) (0/n 2 77- 2 ) + (15/4) ((f>lri 2 7T 2 ) 2 + . . . ] exp (-ti 2 tt 2 I<\>). (2-14) 

H = l 

The repeated occurrence of powers of (/> 1/2 in these relations is particularly worthy of note. 
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From (2-10) 



3. Evaluation of F r (<£) for r>=2 



= Lim 

</r-*0 



+ (1/4)0- (2/3) (l/VV)</> 3 ' 2 + (l/8)c/> 2 + (1/30) (l/vV)^ 2 + (l/96)(^ 3 + . . . .1 
= Lim y ( i/ M (i- e -^). 



(3-1) 



From the relation [1, p. 581, eq (1) with v = 0] 



JoiPnz) 



y z, and 

■/l(ffn*) _ 



= 1/2, 



by multiplying both sides of this equation by z, and integrating with respect to z from zero to x, the 
expression 






is obtained. From this 



|,K7k)/>^ =(1,4, /> 

and carrying through the indicated integration 

^ JoiPnX) 

From this at x = 



2$Ij&=whi- x *). 



(3-2) 



G 2 (0) = V — 7 — =1/8- 
Multiplying (3-2) through by x and integrating from zero to z, it follows that 



(3-3) 



At z = 1 , then 



cr»<»=£(l//J n ) 4 =.F 2 (0)=l/32. 



(3-4) 



(3-5) 



Combining this result with (3-1), it follows that for small c/> 





FM) 






- + (1/32)- 


-(1/4)</> + (2/3)(1/vV)</> 3 ' 2 - 
-(l/96)4> 3 + .... 


-(l/8)(/> 2 - 


-(1/30)(1/V^)4> 5 ' 2 



(3-6) 



Further examination of this method of successive integrations will show that it can lead to 
the evaluation of all sums of the form F r (0) and G r (0), and 
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Lim V (II pf) exp (-#$) = lim F r (</>). 

4>-*o ^ 0—0 



In particular, application of the technique outlined yields for small $ values 



F 3 (<f>) 

= + (l/192)-(l/32) ( /) + (l/8)(/)--(4/15)(l/V / ^(/) 5 / 2 + (l/24)(/) 3 

+ (1/105X1/ Vtt^ 7/2 + (l/384)</> 4 + . . . . 

(3-7) 

These results do not appear to be well known except for F,(0) = cr ( r) which can be found in 
general by a procedure due to Lord Rayleigh [9]. These quantities are all rational fractions, the 
first few of which are as follows: 

r /r r (0) = o- (r) 

1 1/4 

2 1/32 

3 1/192 

4 11/12,288 

5 19/122,880 

For integer values of r, cr ( r) is half the coefficient of z 2r ~ l in the expansion of J\{z)/Jo{z). In fact, 
this follows by integrating Ji{z)l[z 2r Joiz)] around the contour given in section 2 above. 3 See also, 
the paper by Speigel [10]. 

4. Evaluation of Gi(</>) and G>(0) 

The temperature v(R, t) in an infinite circular cylinder of radius a and thermal diffusivity k, 
as a function of time, t, for the boundary conditions 

v(R,0) = v(a,t)=l 
is given by [7, p. 328] 

v (/?, t) = 1 - 2 f J< f£{ a) exp [- K(fiJaM 

The function has the Laplace transform 



(4-1) 



t iwj? *vi— i --1 hiRvpR) 
L p [v(K, t)\- + p Tr=rl ' 

lo(aVp/K) 

Using the asymptotic expansion for 7 (z), eq (2-6) above, it readily follows that 

Lim L p [v(0, t)] 



(4-2) 



r->0 



1/2 



= m P" 3/4 [l~ (1/8) (V^/a)p~ 1/2 - (7/128) (V~KlaYp- l + . . .] exp (- ap 1 ' 2 ! VJc) (4-3) 
Each of the terms in (4—3) is of the form 

p -(2A-+D/4 exp {—apUt/yfe) (4-4) 

•'This technique for finding of was suggested by F. W. J. Olver of the NBS, in a letter dated March 29, 1962. 
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where k =1,2,3, .... The inverse Laplace transform of this expression is 

2 -i/4+(i/2)& 7r -i/2 r 3/4+ WV*D(H2)-k(a 2 l2Ktyl 2 exp (- a 2 /SKt) (4-5) 

where the D„(z)'s are parabolic cylinder functions [11, 12]. Inserting this expression and 

</> = Kt/a 2 (4-6) 

into (4—3), it follows that for small values 

Lim t>(0, </>) 

- (8/0) 1 / 4 [+ i>_ 1/2 (l/V2^) 

-(1/4V2)4>*£>_3/ 2 (1/V2^j 
-(7/64 </> £>_ 5 /2(l/V2^)+...]exp(-l/8<A). (4-7) 

Carsten and McKerrow [13] have obtained the leading term in (4-7) in a somewhat different form. 
It is known that [11, p. 119] 4 



D- 1 i 2 (x) = (xl2iryl 2 K 1 i 4 (x*/4) 



(4-8) 



where K v (x) is the Bessel function of the second kind with a purely imaginary argument [1, p. 77]. 
From this, then, for small <j> 



iA0, < /.) = (l/V^)e-" 8 *A: 1 /4(l/8^)[l + . . . .] 



(4-9) 



For some purposes this may be an easier form to use than (4-7) — the £1/4(2) values are tabulated 
[13] for < 2(0.1) < 5.0 and also for 2 = 6, 8, and 10. 

Higher terms in Lim iA0, </>) are of interest and to obtain these, use is made of the asymptotic re- 



lation [14, p. 92] as 2-» 00 



DM 



~(z v e-* i l*)[\-v{v-\)(ll2z i ) + (\l2)v{v-\)(v-2)(v-Z){ll2z i y + ...} 
and from this, it follows that for small </> 

£U /2 (1/V2£> 
«(2</>) +W4 e-" 8 *[l-(l/4)A:(A:+2)0 + (l/32)yc(A:+2)(A: + 4)(A; +6)<£ 2 + . . .] 



and 



Now since (cf (4-1)) 



it follows that for small $ 



v(0, <f>) - 2e~ l l**(l - </> + 4</> 2 + . . . . ) 



e-" 5* = (1/2X1 -««),«£)] 



A/SJiOiO 



(1/2) - (1 - 4> + 4</> 2 + . . . ) exp (- l/4</>). 



(4-10) 



(4-11) 



(4-12) 



(4-13) 



4 The expression given in (11, p. 119] is ambiguously, if not incorrectly, printed. Equation (4-8) is the proper form, as can be shown by using the relations [12] 
D-, n (x) = 2-^e-^ )I1 U(li; M\ x*/2) and [12, p. 13, eq 1.8.7] U(Ji; %; xV2) = (xV2ir 2 )'/«e + < 1 /«»* , K ±I/ 4(x«/4). 
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Since 



Jo lW „5 *Vl<W 



integrating (4-13) and using (3-3), 



1 



where 



exp (-/3 2 0) = Lim 



+ (1/8) -(1/2)0 
(1/6)0(7 - 40 + 80 2 + . . .) exp (- 1/40) + (7/24)£i(- 1/40) 
Ei(- X )=- \ {z- l e~ z )dz 



(4-14) 
(4-15) 



is the exponential integral function, described and tabulated in Jahnke-Emde [8]. Asymptotically, 
for large values of the independent variable 



Ei(-x)~-x- 1 e-*[l-(llx) + 2(llx) 2 -6(llxf + ...] 
and introducing this in (4-14), for small 0's 



(4-16) 



C 2 (0) 
« + (1/8) - (1/2)0 + 40 2 ( 1-90 + ...) exp (-1/40) 



(4-17) 



Successive integration in the same way, will yield G r (0) for any r^/3. Thus, for example, 




(4-18) 



Differentiation of (4-13) with respect to leads to 



G o (0) 
■■ 4- (1/40) 2 (1 - + . . .) exp (- 1/40) 



(4-19) 



and further differentiation gives expressions for G- r (0) in the limit of small 0. 



5. Computation of <t (3/2) and 



00 



(5/2) 



A general method for finding 



o-o"-» = F r (0) = J(l/)3 n ) 2 



(5-1) 



where r is an integer ^ 1 has been developed above. In the next section, it is necessary to use 
<To (3/2) and cr (5/2) — i.e., odd rather than even powers of /3 n occur in the sums. In contrast to the 
previously considered case, it does not appear possible to express the sum of the odd powers of fi n 
in closed form. 5 The numerical computation of o" (r/2) is simplified by noting that (1-2) yields very 
good approximations to /3 n " r . Thus, for example, the difference between (l//3 n ) 3 and (l//3 n *) 3 is very 
small when only the first three terms of the asymptotic series are taken, e.g., 



Note the similarity to the sum V l/n r = TT r lv{r), where v{r) is an integer for even values of r; but the sum cannot be reduced to a closed form for odd r [6, p. 140]. 
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n 


d/iSn) 3 


d/j3n*) 3 


1 
2 
3 


0.07190 33694 
.00594 51708 
.00154 30866 


0.07172 48723 
.00594 51433 
.00154 30863 



This suggests writing [15] 



(To 



(r/2) 



- £ (l/j8 n ) r = £ \aipn) r -(llPn*) r ]+^(Wn*Y 

n=l n=l L J /i=i 



(5-2) 



Using the asymptotic relation (1-2), terminated at the 0:3 term, for /3 W *, the sums can be evaluated in 
terms of tabulated functions — the poly gamma function [16, p. 9], the sums of negative powers of the 
integers [16, pp. 244, 304], and Euler numbers [16, p. 276, etc.]. Combining these relations, it is 
found that 



o- (3/2) = 0.0808 8147 



(5-3) 



with an uncertainty of no more than one in the last place given. Similarly, using the same tech- 
niques, it has been found that 



o- <5/2> = 0.0126 5566 



(5-4) 



6. The Sums M r (<f>) and N r (4>) 



The functions F r (4>) and G r (</>) are closely related to two other functions 

M r (</,) - J (1//3,,) 2 -- 1 erfc ((3„V$) 



(6-1) 



^erM^V^) 



(6-2) 



which also occur in some cylindrical geometry problems. Here, the complementary error function 

[17] f „ 

erfc (z) - (2/Vw) exp (-£*)<# (6-3) 



From eq (2-10) for small (/> 



r 4> h 

2 ( W) 

n=l Jo 



exp (—p n 2 z 2 )dz 



r<t>> r i 

+ (l/4)-(l/V^)z + 7Z + — -=z* + — z 4 +. . . 
Jo L 4 12V^ 32 



dz. 



(6-4) 



By the change of variable r) = /3nZ, the expression on the left of this equation can be evaluated, and 
it is found that for small (f> values 
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M 2 (<f>) 
cr (3/2) -(l/2) {l/Vn)4> 1 l 2 + (llir)4>--(ll6) (l/V^)^ 2 
-(1/24) (1/tt)(/> 2 - (1/80) (l/V^)(/) 5 / 2 + . . .. 



SothatM 2 (0) = o- (3/2) . 

Similarly, using (3-6) for small 



M 3 (</>) 

^o- (5/2) -(l/16) (l/V^)(/) 1 / 2 + (l/6) (1/V^)(/)^ 2 

-(1/3) (1/tt)(/) 2 + (1/20) (1/Vtt)0 5/2 + (1/9O) (1/tt)</> 3 
+ (1/336) (l/V^)0 7 / 2 + . . . . 



(6-5) 



(6-6) 



It is clear that this procedure can be used to obtain M r (4>) from the value of F r -i(<f>) in each case. 
Note that, using the asymptotic series [18, p. 129] for erfc (j8 n V0) for large (/> values 



M r (4>) 




- -W" 2 t(W fl — (1/2) (/ft,**)- 1 + (3/4) (W)" 2 

V7T £S L 




-(15/8)03„ 2 </>)- 3 +. . . 


exp(- 


-li,?4>) 


-^fe*-" 


F,#)-(l/2<^)F r+1 (</.) + (3/4^)F,. +2 (^) 




-(15/80 s )F r+3 (0)+. . .. 







(6-7A) 



(6-7B) 



The /V,-(</>) functions are obtained in a similar way, by integration, of the corresponding G r (<f>). By 
this procedure, it is found that for small 4> 




(6-8) 



n-m- 



— ^= (3 - 4</))0" 2 — l -j= </> 7 ' 2 exp (- 1/40) + . 

12 V77 V7T 



(6-9) 
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The leading terms in these expressions can be computed by rapidly convergent series to be 

iVi(0) = 0.26659. . ., Ni(0) = 0.054 972. . .. 
For large 4> values 



NA,<t>) 




l 

V7T n =i 


1 - (1/2) G8» 2 4>)-' + (3/4) (i8„ 2 (/))- 2 - (15/8) (£„ 2 </>)- 3 + . . 


•] 


w exp(-/3n 2 <^) 

j8n 2r+ Vl(i8n) 




=^*- w 


G r+1 (</>) -(l/2<£)G r+2 ((/>) + (3/4<^ 2 )G,. +3 ((A) - (15/84> 3 )G r+4 W>) + . 







(6-10A) 

(6-10B) 
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